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Abstract 

Results are presented on the photoproduction of isolated prompt photons, inclusively and 
associated with jets, in the jp center of mass energy range 142 < W < 266 GeV. The 
cross sections are measured for the transverse momentum range of the photons 5 < Ej, < 
10 GeV and for associated jets with E^ 1 > 4.5 GeV. They are measured differentially as 
a function of E^, E^ , the pseudorapidities rf 1 and rf et and estimators of the momentum 
fractions x 1 and x p of the incident photon and proton carried by the constituents participat- 
ing in the hard process. In order to further investigate the underlying dynamics, the angular 
correlation between the prompt photon and the jet in the transverse plane is studied. Pre- 
dictions by perturbative QCD calculations in next to leading order are about 30% below the 
inclusive prompt photon data after corrections for hadronisation and multiple interactions, 
but are in reasonable agreement with the results for prompt photons associated with jets. 
Comparisons with the predictions of the event generators PYTHIA and HERWIG are also 
presented. 
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1 Introduction 



Prompt photon emission in hadronic interactions is a sensitive probe of QCD dynamics 
and partonic structure, providing complementary information to the study of jet production. 
Although cross sections are smaller in the prompt photon case, an isolated photon at large 
transverse energy £y can be related directly to the partonic event structure. This contrasts with 
jet measurements, where the partonic structure is obscured by the non-perturbative hadroni- 
sation process. Furthermore, experimental uncertainties associated with the transverse energy 
measurement are smaller for the electromagnetic showers initiated by photons than for the mea- 
surement of jets of hadrons. This paper presents results from the HI experiment at HERA on 
the photoproduction of isolated prompt photons, both inclusively and in association with a jet. 
Here photons are called "prompt" if they are coupled to the interacting quarks (see Fig. [T|), in 
contrast to photons which are produced as decay products of hadronic particles. 

Next to leading order (NLO) perturbative QCD (pQCD) calculations have proved broadly 
successful in reproducing measured jet production rates from various colliders, provided hadro- 
nisation corrections are applied. In contrast, discrepancies have been observed between data 
on prompt photon production and NLO pQCD calculations in pp, pp, pN (see e.g. [1]) and 
in 7p [2] interactions. These discrepancies can be reduced by introducing intrinsic transverse 
momentum k T to the incoming partons of the proton [1,3] or by soft gluon resummations [4]. 
NLO pQCD calculations of prompt photon production in photoproduction [5-9] are available, 
in which the incident photon interacts with the partons of the proton either directly (Fig.[lk,b) 
or via its "resolved" partonic structure (Fig.[T]:,d). This paper investigates the extent to which 
fixed order NLO calculations are able to describe the new data. The data are also compared 
with the predictions of the event generators PYTHIA [10] and HERWIG [11] and with data on 
inclusive prompt photon production from the ZEUS collaboration [2]. 

2 Strategy of Prompt Photon Measurement 

The photoproduction process is initiated by quasi-real photons, which are produced in small 
angle ep scattering, where the scattered electron 1 escapes into the beam pipe. The requirement 
of a significant energy loss in the electron beam direction, together with the condition that no 
electrons are found in the selected events, suppresses contributions from neutral current (NC) 
deep inelastic scattering (DIS). 

Photons are identified in the HI liquid argon (LAr) calorimeter [12] as compact electro- 
magnetic clusters with no track associated to them. The main experimental difficulty is the 
separation of the prompt photons from hadronic backgrounds, in particular tt° meson decays, 
since at high energy the decay photons cannot be reconstructed separately in the calorimeter. 
Since the ir° mesons are predominantly produced in jets, this background is strongly reduced 
by requiring the photon candidates to be isolated from other particles. 

After the selection cuts described below, the tt° background is still of similar size to the 
prompt photon signal. To extract the signal, different shower shape variables are combined 

'The term "electron" is used for both electrons and positrons. 
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to form a discriminator function which is fitted with a sum of contributions from simulated 
photons, 7T° and 77 mesons. The fit is done double-differentially in bins of transverse energy Ej, 
and pseudorapidity rf 1 . 2 

The number of selected prompt photon events is corrected for detector effects by detailed 
simulations of prompt photon production in the HI detector, using the event generators PYTHIA 
and HERWIG. 

3 Event Selection 

The data were collected in the years 1996-2000 with the HI detector [13] at HERA in data 
taking periods where electrons or positrons with energy E e = 27.6 GeV collided with protons 
of energies E p = 820 GeV or E p = 920 GeV. The data correspond to an integrated luminosity 
of 105 pb _1 of which 28.8 pb _1 and 61.3 pb _1 were recorded in e + p interactions at center of 
mass energies y/s = 301 GeV and y/s = 319 GeV, respectively, and 14.9 pb _1 were recorded 
in e~p interactions at y/s = 319 GeV. 

The events are triggered by compact energy depositions in the LAr calorimeter, consistent 
with electron or photon showers. The trigger efficiency is ps 60% at photon energies of 5 GeV, 
reaching 100% at w 12 GeV. 

The main requirements for the event selection are the following [14]. 

• A compact electromagnetic energy cluster, consistent with a 7 shower, is reconstructed 
in the LAr calorimeter in the range —1 < rf 1 < 0.9 and E], > 5 GeV. For the data at 
y/s = 301 GeV a threshold Ej, > 7 GeV is required. 3 The if range corresponds to the 
central barrel region of the LAr calorimeter [12]. 

• No track is allowed to point to this cluster within a distance of 25 cm in the plane trans- 
verse to the track at the calorimeter surface. 

• Events with electron candidates in the LAr calorimeter or in the backward calorimeter 
SPACAL [15] are rejected. This restricts the virtuality of the exchanged photon to Q 2 < 1 
GeV 2 and suppresses contributions from radiative DIS and QED Compton processes. 

• At least two tracks are required in the central tracker [13, 16], which covers the angular 
range 20° < 9 < 160°. This cut assures good reconstruction of the event vertex, which is 
required to be within ±35 cm in z of the nominal vertex position. 

• The inelasticity y = W 2 /s, where W is the ^p center of mass energy, is evaluated as 
y = Yl {E — Pz)/2E e . Here the sum runs over all detected final state particles. The 
required range 0.2 < y < 0.7 corresponds to 142 < W < 266 GeV for E p = 920 GeV. 
The cut at high y reduces NC DIS background. The cut at low y removes beam gas 
background. 

2 The pseudorapidity r\ of an object with polar angle 8 is given by 77 = — In tan(#/2), where 9 is measured 
with respect to the z axis given by the proton beam direction. Transverse energies are also measured with respect 
to this axis unless otherwise specified. 

3 The data at y/s = 301 GeV were taken prior to the upgrade of the LAr electronics, after which lower trigger 
thresholds became possible. 
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• The 7 candidate is required to be isolated. The transverse energy, E^? ne , in a cone around 
the 7 candidate, given by distances below 1 unit in the (77 — <p) plane, is required to be 
less than 10% of Ej, . A further cut to remove non-prompt photon background is based 
on the shower shape, as described in section |4j 

• Associated jets are reconstructed using the inclusive kr algorithm [17] with a distance 
parameter D = I'm. the (77 — <p) plane. The jets are selected in the jet energy and pseudo- 
rapidity ranges E 3 ^ > 4.5 GeV and — 1 < rf 6 * < 2.3, respectively. If more than one jet 
is found, only the jet with the highest is considered. 



4 Signal Extraction 

The prompt photon signal is extracted using a shower shape analysis based on the expectation 
that showers initiated by photons are typically narrower than n° or rj initiated showers, with 
less energy deposited in the first calorimetric layer on average. This procedure exploits the fine 
segmentation of the electromagnetic section of the LAr calorimeter, which has transverse cell 
sizes varying between about 5x7 cm 2 and 7 x 13 cm 2 in the central barrel region. This part 
of the calorimeter has three layers in depth, corresponding to « 20 radiation lengths. The first 
layer has a thickness of about 3 radiation lengths. 

Three observables are used to discriminate against background. The mean transverse shower 
radius is given by R = £\ r^,;/ £\ e is where r* is the transverse distance of cell i with energy 
density S{ measured with respect to the axis from the event vertex to the center of gravity of the 
7 candidate cluster. The shower hot core fraction (HCF) is the largest energy fraction of the 
candidate shower which is contained in 4 or 8 contiguous cells (depending on the calorimeter 
granularity) including the cell of highest energy. Finally the first layer fraction (FLF) is the 
energy fraction of the shower contained in the first, i.e. closest to the beam, layer of cells of the 
calorimeter. The observables R and FLF are expected to be smaller and HCF to be larger for 
the prompt photon showers than for background showers. 

To discriminate between photon and background showers, probability densities for the three 
observables R, HCF and FLF are determined by simulation of photons and ir° mesons [14]. 
The products of these three densities are used as likelihood functions. For each measured event 
a discriminator (d) is formed by the likelihood for photons divided by the sum of the likelihoods 
for photons and n° mesons. The discriminator d produces larger values for prompt photons than 
for the pairs of photons from %° or 77 meson decays. This can be seen in Fig.|2^, where the mea- 
sured distribution of d is shown together with the prompt photon and background components 
estimated from the PYTHIA Monte Carlo simulation. The lowest bin in Fig.Ek contains mostly 
background 90%) which is composed of 65% 7r°, 30% 77 and 5% other particles. After a 
cut d > 0.125, the remaining background (~ 50%) is composed of 94% 7r°, 5% 77 and about 
1% other particles. Fig. shows that the distribution in the discriminator d of signal plus 
background is well predicted by PYTHIA. 

The contribution of prompt photons is determined by maximum likelihood fits of simulated 
photon, 7T° and 77 meson discriminator distributions to the data distribution for d > 0.125. Each 
measurement presented in section 6 is obtained by summing the results of such fits performed 
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independently in 6 x 6 bins of rf and E^. In this procedure, only the 77/71-° ratio (on average 
5% after the selection requirements) is taken from PYTHIA. 

The measured distributions of R, HCF and FLF for the full rp and Ej, range are shown 
in Fig.|2j5-d together with the contributions of photons and 7r° and 77 background, the normali- 
sations of which are taken from the fits to the discriminator distributions described above. The 
data distribution is well described by the extracted signal and background components. The 
discrimination between signal and background becomes weaker at high Ej, , where the R and 
HCF distributions of ir° mesons and photons become more similar to each other. Therefore 
events with E^ > 10 GeV are not included in the results presented below. 

5 Monte Carlo Generators and Corrections to the Data 

The event generators PYTHIA [10] and HERWIG [11] are used 4 to correct the observed event 
yields for apparatus effects by means of a full simulation of the HI detector. The average cor- 
rections from the two generators are applied. Both generators are based on leading order (LO) 
QCD matrix elements and leading log parton showers. Hadronisation is provided in PYTHIA 
by Lund string fragmentation [18] and in HERWIG by the decay of colourless parton clusters. 
Both generators model additional soft remnant-remnant interactions, termed multiple parton 
interactions (m.i.) in the following. The GRV(LO) [19] parton densities are used for the pho- 
ton and the proton. In contrast to HERWIG, PYTHIA simulates radiation of photons from the 
electron line and photon production via fragmentation of final state quarks and gluons in di-jet 
events. In PYTHIA the parameter describing the intrinsic kr of initial state partons in the pro- 
ton is k = 1 GeV (default) leading to < k T >= 0.9 GeV. HERWIG predictions are shown for 
the default value k T — 0. 

A correction factor of about 1.04 is applied to the y/s = 301 GeV data in order to combine 
with the data at y/s = 319 GeV. The final results are presented for y/s = 319 GeV. Background 
from DIS events where the scattered electron fakes the prompt photon signatures due to tracker 
inefficiency leads to a subtraction of 3% in the lowest rf bin at high E\ , and is negligible 
otherwise. Following the selection criteria described in section 3, other sources of background 
are negligible. 

6 Systematic Uncertainties 

The following sources of systematic errors are considered. 

• The dominant systematic errors are due to possible imperfections in the simulation of the 
shower shapes. The quality of the simulated distributions of R, HCF and FLF is tested 
by comparing simulations of electrons with electron candidates from NC DIS events. 
Differences between the data and the simulation in these distributions result in errors on 
the prompt photon cross sections ranging from ±10% to ±20%. 

4 PYTHIA 6.15/70 and HERWIG 6.1 are used with default multiple interactions (MSTP(82)=1 for PYTHIA, 
PRSOF=l, BTCLM=1 for HERWIG). 
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• The uncertainties on the calorimeter electromagnetic and hadronic energy scales (0.7% 
to 1.5% and rs 3% respectively) contribute errors of about 5% for the inclusive cross 
sections. For the case with associated jets, the hadronic energy uncertainty gives rise to 
uncertainties of about 10%. 

• The trigger efficiency is determined using independent triggers with an uncertainty of 3%. 

• The model dependence, of the corrections for detector effects is quantified as half the dif- 
ference (< 6% in most cases) between the correction factors obtained with PYTHIA and 
HERWIG. Within these uncertainties the results are insensitive to reasonable variations 
of the Ej, dependence and of the underlying event activity in PYTHIA. 

An overall normalisation uncertainty of ±1.5% on the luminosity measurement is not in- 
cluded in the results. The total systematic errors are obtained by adding the different systematic 
errors in quadrature. For further details see [14]. 



7 Results 

The results are presented in tables 1 to 3 and figures|3]to |5]as bin averaged ep cross sections in 
the kinematic region defined by 

= 319 GeV, 0.2 < y < 0.7, Q 2 < 1 GeV 2 

and 5 < E? T < 10 GeV, -1< r/ 7 < 0.9, E^ ne < 0.1 • E% . 

The inner error bars on the data points in the figures indicate the statistical errors as obtained 
from the shower discriminating fits. The full error bars also contain the systematic errors added 
in quadrature. 

7.1 Inclusive prompt photons 

Differential cross sections da/dE^ and da /drf for inclusive prompt photon production are 
shown in Fig. |3] and are compared with the predictions of the PYTHIA [10] and HERWIG [11] 
event generators. The cross sections are reasonably described in shape, but the predictions by 
PYTHIA (HERWIG) are low by about 40% (50%) in normalisation. Photons from fragmenta- 
tion in di-jet events are only treated in the PYTHIA calculation, which explains the difference 
from HERWIG. The figure shows the full PYTHIA prediction and separately the contribution 
of direct photon interactions only. The PYTHIA prediction is also shown without multiple in- 
teractions. The predictions at < r[ 1 < 0.9 are about 25% higher without multiple interactions, 
showing that the cross section is reduced by the soft underlying event activity, as expected [6] 
due to the isolation cone condition E^ ne = 0.1 • E^. Fig.|3]also shows a comparison with the 
results of the ZEUS collaboration [2]. 5 The two measurements are consistent. 

5 The ZEUS data are obtained in the somewhat different kinematic region 0.2 < y < 0.9, —0.7 < rf 1 < 0.9, 
•s/s = 301 GeV and are adjusted to correspond to the HI conditions using the NLO calculation [6]. 
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The results are further compared in Figs.|3^,b with the NLO pQCD calculations by Fontan- 
naz, Guillet and Heinrich (FGH) [6] and Krawczyk and Zembrzuski (K&Z) [8,9]. The two 
calculations are similar, the main difference being that only FGH apply higher order corrections 
to the resolved photon processes (Fig. [lt,d). The final state photon can be emitted in a hard 
partonic process or can be produced in the fragmentation of a quark or gluon. The photon isola- 
tion requirement in the cross section definition suppresses the latter contribution considerably. 
According to the calculations, the contribution to the cross sections from fragmentation is typi- 
cally 10%. Both calculations use the photon and proton parton density functions AFG [20] and 
MRST2 [21], respectively, and BFG [22] fragmentation functions. The transverse energy Ej, is 
used for the renormalisation and factorisation scales. In order to obtain a realistic comparison of 
data and theory, the NLO calculations are shown after correction from the parton to the hadron 
level including the effects of multiple interactions (m.i.). These corrections are obtained using 
the average of PYTHIA and HERWIG, taking half the difference as an uncertainty estimate. 
The resulting uncertainty is typically 3% as shown by the outer error bands in Figs. |4ji,b. 

The FGH (K&Z) NLO calculations are typically 30% (40%) below the data in most of the 
Ej, and rp ranges presented in Figs. Ek,b, if the corrections for hadronisation and m.i. are 
applied. The comparison with the parton level result of FGH shows that the correction factors 
are largest at high rp ', where resolved photon interactions contribute most. Only here (r/ 7 > 0.6) 
the corrections improve the agreement with the data. The NLO corrections are substantial, with 
the NLO/LO ratio increasing from 1.2 to 1.4 for FGH with increasing rp . The effect on the NLO 
calculations of scale variations is rather small as shown by the inner error bands in Fig. 0K,b, 
which indicates the effect of the variation of the renormalisation and factorisation scales from 
0.5 • E% to 2 • E%. The predictions are about 10% larger on average if the GRV parton density 
and fragmentation functions [19,23] are used. 

7.2 Prompt photons with jets 

Cross sections for the production of a prompt photon associated with a jet are presented in 
Figs.Ht,d as a function of the variables Ej, and rp and in Fig.|5]as function of E^, rf et ; , x^° and 

x^°. The estimators x^° = E^(e^ et + e~^)/2yE e and x L ° = E^(e^ et + e^)/2E p are taken 
for the momentum fractions of constituents of the incident photon and proton, respectively, 
participating in the hard process. These observables make explicit use only of the photon energy, 
which is better measured than the jet energy. They are thus most easily interpreted in the 
leading order (LO) approximation where the outgoing partons from the hard interaction, and 
correspondingly the reconstructed photon and jet, have equal transverse momenta. The use of 
the variable x^° was recommended in [7,20] to reduce infrared sensitivity. The variable x^° is 
discussed e.g. in [24]. 

The data are compared with the pQCD calculations. In contrast to the inclusive case, both 
NLO calculations are consistent with the data in most bins. The NLO/LO correction ratios for 
dcr/d?7 7 are more moderate than in the inclusive case, ranging from about 0.95 to 1.25 with 
increasing rp for FGH (Fig.UJi), but are still large in some bins of other distributions, as shown 
in Fig. |5J The hadronic and m.i. corrections, which are applied for both NLO calculations, 
improve the description of the data only in some regions, such as rp > 0.6, if et < —0.3 and 
x^° < 0.6. The x^° distribution (Fig. Efc) is particularly sensitive to the photon structure 
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function. Using the GRV parameterisation [19,23], the K&Z prediction increases by about 
4% for > 0.85 and by about 20% at < 0.85 where resolved photon contributions 
dominate, leading to a somewhat improved description of the data. 

Further understanding of the dynamics of the process and in particular of the effect of higher 
order gluon emissions beyond the diagrams shown in Fig.UJmay be obtained from the transverse 
correlation between the prompt photon and the jet. The distribution of the component of the 
prompt photon's momentum perpendicular to the jet direction in the transverse plane, 

px_ = | ft? x Vt* I / I Vt 1 I = ' sin(A0) , 

where A0 is the difference in azimuth between the photon and the jet, is determined by higher 
order effects as p± is zero at leading order, where the prompt photon and the jet are back-to-back 
in the transverse plane. 

The normalised pj_ distribution is shown in Fig. |6] separately for the regions > 0.85 and 
x^° < 0.85, where direct and resolved photon induced processes dominate, respectively (con- 
tributing about 80% in each case according to the FGH calculation). The FGH NLO prediction 
gives a better description of the p±_ distributions than K&Z. This is due, at least in part, to the 
differences between the QCD corrections for resolved photon interactions in the two calcula- 
tions. For x^° > 0.85, the data are quite well described by PYTHIA, whereas the HERWIG 
prediction is somewhat harder than that of the data. For x^° < 0.85, the measured distribution 
is reasonably well described by both Monte Carlo models, except in the region of p± ~ 5 GeV. 
The p± distribution at large x 7 has previously been used [3] to extract information on an in- 
trinsic k t of the initial state partons of the proton. However, the large differences between the 
predictions of the various NLO calculations and Monte Carlo models in the present comparisons 
do not allow a reliable statement to be made. 



8 Conclusions 

The photoproduction of prompt photons, both inclusively and associated with jets, is studied. 
The measured rf and Ej, distributions of the inclusive prompt photons are reasonably well 
described in shape by NLO pQCD calculations, but after corrections for hadronisation and 
multiple interactions the predictions are 30%— 40% below the data. The cross sections predicted 
by the PYTHIA and HERWIG event generators describe the data distributions well in shape 
with normalisations that are low by about 40% — 50%. 

For prompt photons associated with a jet, the data are somewhat better described by the NLO 
calculations including corrections for hadronisation and multiple interactions. This, together 
with the fact that the NLO corrections are smaller on average than in the inclusive case, suggests 
that contributions beyond NLO are less important if an energetic jet is selected together with 
the prompt photon. 

The distribution of p±, the component of the prompt photon's momentum perpendicular to 
the jet direction in the transverse plane, is sensitive to effects beyond LO. The PYTHIA gener- 
ator describes the normalised p± distributions quite well, whereas HERWIG predicts too hard 
a p± distribution at large x^° , where direct photon interactions dominate. Particularly at low 
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, the p\_ distribution is better described by the NLO calculations if NLO QCD corrections 
are also applied in the case of the resolved photon interactions. 
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E% (GeV) 


&r/dE} (pb/GeV) 


h.c.+ m.i 


5.0-5.8 


17.0 ±1.8± 2.9 


0.75 


5.8-6.7 


11.6 ±1.2± 1.9 


0.78 


6.7-7.5 


8.3 ±0.9± 1.3 


0.84 


7.5 - 8.3 


6.67 ±0.69±0.85 


0.88 


8.3 - 9.2 


5.46 ±0.55±0.93 


0.91 


9.2 - 10.0 


2.42 ±0.48± 0.57 


0.91 


if 1 


da/drf (pb) 


h.c.+ m.i. 


-1.0 - -0.7 


26.8 ±2.6±4.2 


0.94 


-0.7 - -0.4 


29.1 ±3.4± 5.2 


0.87 


-0.4 - 0.0 


26.9 ±3.0± 4.3 


0.84 


0.0 - 0.3 


24.7 ±2.7± 3.4 


0.75 


0.3 - 0.6 


21.2 ±1.9± 3.5 


0.71 


0.6 - 0.9 


11.0 ±1.6± 2.2 


0.65 



Table 1: Inclusive prompt photon differential cross sections da /dE^ for — 1 < rp < 0.9 and 
da/drf for 5 < E% < 10 GeV with = 319 GeV and 0.2 < y < 0.7. The first error is 
statistical, the second systematic. The correction factors applied to the NLO calculations for 
hadronisation and multiple interactions (h.c. + m.i.) are also given. 
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El (GeV) 


da/dEl (pb/GeV) 


h.c.+ m.i 


5.0-5.8 


8.9 ±1.2± 2.4 


0.77 


5.8 - 6.7 


6.24 ±0.90± 1.02 


0.80 


6.7 - 7.5 


6.10 ±0.77± 0.97 


0.85 


7.5 - 8.3 


5.28 ±0.64± 0.71 


0.89 


8.3 - 9.2 


3.98 ±0.51±0.71 


0.93 


9.2 - 10.0 


2.30 ±0.48± 0.52 


0.91 


rf 


da/drf (pb) 


h.c.+ m.i. 


-1.0- -0.7 


16.2 ±2.1± 2.8 


0.95 


-0.7 - -0.4 


18.2 ±2.5± 3.2 


0.89 


-0.4 — 0.0 


18.2 ±2.4± 2.9 


0.87 


0.0 - 0.3 


13.7 ±2.0± 2.2 


0.77 


0.3 - 0.6 


14.5 ±1.5± 2.6 


0.76 


A f. A A 

u.d — u.y 


A C Q_l_ 1 C. 
O.O ±l.d± 1.0 


U. / 1 


(GeV) 


do-/dE^ (pb/GeV) 


h.c.+ m.i. 


4.5 - 6.7 


5.97 ±0.54± 0.91 


0.82 


6.7 - 8.8 


3.92 ±0.42± 0.72 


0.84 


8.8 - 11.0 


2.07 ±0.26± 0.44 


0.83 




da/d^ et (pb) 


h.c.+ m.i. 


-1.0- -0.3 


7.5 ±1.2± 1.4 


0.77 


-0.3 - 0.3 


12.4 ±1.2± 1.6 


0.88 


0.3 - 1.0 


12.3 ±1.3± 2.0 


0.86 


1.0 - 1.6 


6.58 ±0.87± 1.10 


0.78 


1.6 - 2.3 


2.81 zbl.Olzb 0.94 


0.91 


m 

, T LU 

•s 


da/dx^° (pb) 


h.c.+ m.i. 


0.1 - 0.3 


7.8 ±2.2± 2.5 


0.50 


0.3 - 0.6 


16.4 ±2.8± 2.6 


0.65 


0.6 - 0.9 


39.9 ±4.1± 7.4 


1.1 


0.9 - 1.1 


49.3 ±4.1± 7.4 


0.88 


x p 


da/dx^ u (pb) 


h.c.+ m.i. 


0.0018 -0.0034 


247 ±149± 51 


0.84 


0.0034 - 0.0063 


2350 ±300± 350 


0.80 


0.0063 - 0.0120 


1900 ±180± 290 


0.86 


0.0120 - 0.0220 


640 ±70± 113 


0.81 


0.0220 - 0.0400 


167 ±31± 39 


0.84 



Table 2: Prompt photon cross sections with an additional jet requirement (Ej? > 4.5 GeV, 
-1 < rf et < 2.3) differential in E%, rf, E 3 T et , r] jet , and for -1 < rf < 0.9 and 
5 < E^ < 10 GeV with s/s = 319 GeV and 0.2 < y < 0.7. The first error is statistical, the 
second systematic. The correction factors applied to the NLO calculations for hadronisation 
and multiple interactions (h.c. + m.i.) are also given. 
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P± (GeV) 


l/o- da/dj9 ± (GeV~ 1 ), < 0.85 


0-2 
9 — 4 

4-6 
6-8 


0.216 ±0.030± 0.015 

0.124 ±0.019± 0.011 
0.0225 ±0.0081± 0.0077 


P± (GeV) 


l/o- d ( j/dp ± (GeV- 1 ), > 0.85 


0-2 
2-4 
4-6 


0.420 ±0.033± 0.024 
0.061 ±0.017± 0.014 
0.0054 ±0.0078± 0.0026 



Table 3: Normalised cross sections differential in the prompt photon momentum component 
perpendicular to the jet direction in the transverse plane, for < 0.85 and > 0.85 with 
5 < < 10 GeV, -1 < jf < 0.9, E J T et > 4.5 GeV and -1 < rf et < 2.3 with yfi = 319 GeV 
and 0.2 < y < 0.7. The first error is statistical, the second systematic. 



a) 



c) 




Figure 1 : Examples of leading order diagrams producing prompt photons. Direct photon inter- 
actions a), b) and resolved photon interactions c),d). 
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Figure 2: a) Distribution in the discriminator d, with the analysis cut indicated, for data (solid 
points), and PYTHIA normalised to the data (solid line), with prompt photons (dashed) and 
the sum of the background contributions (dotted), b) to d) Distributions of the mean transverse 
shower radius R, the hot core fraction HCF and the first layer fraction FLF for all selected 
photon candidates (data points). The contributions determined for photons (dashed lines), back- 
ground (7T° + r], dotted lines) and the sum (solid lines) are also shown. 
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Figure 3: Inclusive prompt photon differential cross sections dcr/dE^ for — 1 < rp < 0.9 (a) 
and da/dr/ 7 for 5 < < 10 GeV (b) with y/Ts = 319 GeV and 0.2 < y < 0.7 compared with 
the predictions of HERWIG (dashed line) and PYTHIA including multiple interactions (full 
line). The contribution of direct interactions is shown separately (dashed-dotted line). The full 
PYTHIA prediction without multiple interactions (dotted line) is also shown. ZEUS data [2] 
are shown adjusted to correspond to s/s = 319 GeV, 0.2 < y < 0.7 and — 1 < rp < 0.9. 
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Figure 4: Inclusive prompt photon cross sections (a,b) for — 1 < rp < 0.9 and 5 < E^, < 
10 GeV with y/s = 319 GeV and 0.2 < y < 0.7, and with an additional jet requirement 
(E J T et > 4.5 GeV, -1 < rj jet < 2.3) (c,d). The data are compared with NLO pQCD calculations 
(K&Z [9], dotted line, and FGH [6], solid line). The NLO results are corrected for hadronisation 
and multiple interaction (h.c. + m.i.) effects (see text). The outer error bands show the estimated 
uncertainties on these corrections for the example of FGH, added linearly to the effect of a 
variation of the renormalisation and factorisation scales in the NLO calculation from 0.5 • 
to 2 • Ej, (inner band). The FGH results are also shown without corrections for h.c. and m.i. at 
NLO (dashed) and LO (dashed-dotted). 
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Figure 5: Prompt photon cross sections with an additional jet requirement (E^ > 4.5 GeV, 
— 1 < rf et < 2.3) differential in E^, rf et , and x^°. The data are compared with NLO 
pQCD calculations (K&Z [9], dotted line and FGH [6,7] solid line). The NLO results are 
corrected for hadronisation and multiple interaction (h.c.+ m.i.) effects (see text). For the error 
bands see Fig.|4j The FGH results are also shown without corrections for h.c. and m.i. at NLO 
(dashed) and LO (dashed-dotted). 
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Figure 6: Distribution of the prompt photon momentum component, perpendicular to the jet 
direction in the transverse plane, for < 0.85, a) and c), and > 0.85, b) and d). In a) 
and b) the data are compared with PYTHIA (solid line) and HERWIG (dashed line). In c) and 
d) the data are compared with NLO pQCD calculations (K&Z [9], dotted line, and FGH [6,7], 
solid line). The NLO results are corrected for hadronisation and multiple interactions. For the 
error bands see Fig.|U 
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